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H ' We use a continuous-time path integral to obtain the semiclassical prop- 

agator for minimal-spread spin coherent states. We pay particular attention 

to the "extra phase" discovered by Solari and Kochetov, and show that this 

correction is related to an anomaly in the fluctuation determinant. We show 

that, once this extra factor is included, the semiclassical propagator has the 

correct short time behaviour to O(T^), and demonstrate its consistency under 

dissection of the path. 
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I. INTRODUCTION 



Coherent-state path integrals for spin were introduced by Klauder |Q, and by Kuratsuji 
and Suzuki 0. Related phase space path integrals were introduced by Jevicki and Papani- 
colaou 0] , and by Nielsen and Rohrlich 0] . For a review see . These path integrals have 
attracted attention in connection with geometric quantization [0], and for providing exam- 
ples hinting at possible infinite-dimensional extensions of the Duistermat-Heckman theorem 
1^ on conditions for the exactness of the stationary phase approximation @JT^. Perhaps 
their most significant practical applications, however, have been in computations of spin 
tunnelling in the semiclassical limit. Here the spin path- integral formalism gives a good 
qualitative description of the tunneling process [p!lH T3|, including the simplest and most 



vivid picture of the topological quenching of spin tunneling []I4| that has recently been seen 
in the magnetic molecule Feg [l^. When we require precise quantitative results, however, 
the spin coherent-state path integral runs into problems: A straight forward application of 
instanton methods to compute the tunnel splitting yields answers that are incorrect 

beyond the leading exponential order |18[. A full derivation of the splitting, including the 
correct prefactor, has only recently been provided by Belinicher, Providencia and Providen- 
cia |T^ . These authors showed that the continuum limit of the discrete path integral is rather 
delicate, and in their computation the simplicity of the instanton method is lost. These dif- 
ficulties have lead to the spin path integral acquiring a reputation for being unreliable — or. 



even worse, being meaningful only in its discrete-time form pO[. Many workers in the field 



have sought alternatives to path integrals such as discrete WKB methods pl|-P3| . 

This paper is intended to effect a rehabilitation of the continuous-time spin coherent- 
state path integral. We advertise and explain the origin of a previously discovered, but 
largely unknown, correction to the naive form of the semiclassical propagator. This "extra 



phase" was obtained by Solari ||2J] as a result of a careful evaluation of the discrete path 
form of the path integral. It also appears, as a product of a manipulation, apparently 
carried out for convenience, in a paper by Kochetov ||25|. We derive it here by pointing out 
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that the functional determinant resulting from the fluctuation integral about the classical 
path possesses an anomaly. Regulating the determinant in a manner consistent with the 
underlying causal structure leads to the extra contribution. 

The structure of the paper is as follows: In section two we review spin coherent states 
built on highest- or lowest-weight spin-j states. We focus primarily on their holomorphic 
properties. In section three we review the properties of the classical action that appears 
in the path integral for spin, stressing the importance of boundary terms in avoiding the 
over-determination problem. In section four we compute the gaussian integral over small 
fluctuations about the classical path, and obtain the extra-phase correction. In section 
five we verify that, once the extra contribution is taken into account, the semi classical 
propagator has the correct short time behaviour. This verification is immediate at first order 
in T, but the agreement between our expression and the exact result at O(T^) provides a 
significant test of the correctness of our result. In section six we check the consistency of the 
expression for the propagator under the dissection of the path. We find that our semiclassical 
propagator does not pass this test unless we repartition terms between the exponent and 
the prefactor. This forces us to regard the large parameter in the semiclassical expansion 
as being j + 1/2, rather than j. As a byproduct, this observation resolves the mystery 
of the divergent normalization factor that appears in most treatments of the semiclassical 
propagator. Finally, in section seven, we compute the semiclassical propagator for the 
hamiltonian H = z/J|. We confirm that our expression obtains the correct leading and next 
to leading terms in the large-j expansion. 



II. SPIN COHERENT STATES 



We define a family of spin coherent states by 



\z) = exp{z.J+)\j,-j). (2.1) 

These states are not normalized, but have the advantage of being holomorphic in the pa- 
rameter z. Consequently, matrix elements such as {z'\0\z) will be holomorphic functions of 



the variable and anti-holomorphic functions of the variable z' . 



The inner product of two of these states is 



(2.2) 



and the left eigenstates (j, m| of and J^, have coherent-state wavefunctions 



= {3,m\z) = J- -— — — -z^ 

\ [J — my. [J + my. 



(2.3) 



This means that a general element of the spin-j Hilbert space may be represented by a 
polynomial in z of degree n <2j. 

As with any family of generalized coherent states derived from a unitary irreducible 
representation of a compact group, Shur's lemma provides us with an overcompleteness 
relation. In the present case this reads 



Here 2j + 1 appears because it is the dimension of the representation. The symbol (Pz is 
shorthand for dxdy, and the factor 1/(1 + zzy combines with this to make the invariant 
measure on the coset SU{2)/U{1). This coset is, of course, the two-sphere, 5"^, equipped 
with stereographic coordinates. The south pole, corresponding to spin down, is at 2; = 0, 
while the north pole, spin up, is at 2; = 00 — the one-point compactification of the complex 
plane. The remaining factor in the measure, 1/(1 + zzy^, serves to normalize the states. 

The wavefunctions ipl^^z) are singular at the north pole, z — 00. Indeed there is no 
actual state \oo) because the phase of this putative limiting state would depend on the 
direction from which we approach the point at infinity. We may, however, define a second 
family of states 




(2.4) 



z)^ = exp(^J_)|j, j). 



(2.5) 



and form the wavefunctions 




(2.6) 
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These states and wavefunctions are well defined in the vicinity of the north pole, but singular 
near the south pole. 

To find the relation between ■ip^'^^z) and ijj^^\z) we note that the matrix identity 

(2.7) 

coupled with the faithfulness of the spin-| representation of SU{2), implies the relation 

exp{zJ+)w = exp exp (2.8) 

where w = exp(i7rJ2) is the generator of the Weyl group of SU{2). We also note that 



"1 


z 




" 1" 




1 0" 




'-Z " 




"1 







1_ 




.-1 o_ 




1. 




. -z'^ . 







1 



w\3j) = (-I)^^Ij, -j), -j) = \ 



(2.9) 



Thus, 



{j,m\e''^+\j, -j) 

i-iy^ {j,m\e'^+w\j,j) 

(-lf^j,m\e^-'J-{-zf^e-^-''^\j,j) 



(2.10) 



The coherent-state wavefunctions V'm'* i^m "^^Y therefore be regarded as composing a 
single global section, ipm-, of a holomorphic line bundle with transition function z^^ relating its 
components ip^i^z), and V'mHC = 1/-2) in the two coordinate patches. It is the requirement 
that the transition function and its inverse be holomorphic and single valued in the overlap 
of the coordinate patches that forces 2j to be an integer. In the sequel, all coherent states, 
unless otherwise specified, will drawn from the first family, \z). 

The above construction is an example of the Borel-Weil realization of representations 



of compact groups as sections of holomorphic bundles It serves as the paradigm for 

the more general theory of geometric quantization |[7|,|25|]. Because global analyticity is 
characteristic of the minimal-spread coherent states built on highest- (or lowest-) weight 



states, and also serves (via the transition function) to specify the Hilbert space, it is a 
property that should be maintained order-by-order in any approximation scheme. 

For physical interpretations we must normalize the coherent states. This we do by 
multiplying them by 

N{z,z) = (l + zz)-^. (2.11) 

For example, 

N'{z\Js\z)=j^^, and N'{z\J^\z) = J^. (2.12) 
zz + 1 zz + 1 

If we recall the connection between stereographic and spherical polar coordinates, 

z = e-'^cot^, (2.13) 

we see that 

j^^^=jcos^, and = je-'^sine. (2.14) 
zz + 1 zz + 1 

We also note that 

1^ /zz-l\^ 1 



j(j-^)cos2^ + |. (2.15) 



Similarly p| 



N^z\Jl\z) = j{j - ^) sin2 e cos2 + |, 

N^z\J^\z) = j{j - ^) sin^ e sin^ + |. (2.16) 

Thus N'^{z\P\z) = j{j + 1), as it should. 

The normahzed wavefunctions N{z, z)^^\z) have their maximum amplitude on the lines 
of latitude 

kr = k^r = ^ (2.17) 

J — m 



corresponding to the polar angle 9m = cos ^ m/j. Note that 



P - 1 

N'izMlzm) = j 'r"' , , = m. (2.18) 



z P + 1 



The variance, in terms of m, is given by 



N'{J!) - NVsf) = \j (l - (^)') = - cos^ 0). (2.19) 

Since m ^ j cos 6', the normalized wavefunctions have zonal spread ~ 1/ ^/]. As j 
becomes large the quantum spin becomes more localized, and more classical. 



III. SPIN ACTION 

We wish to find a semiclassical approximation for the propagator 

^(C/,0,T) = (C/|e-^^^|0) (3.1) 

in the form 

A'scKC/, Cu T) = A'reduced " GXp {^c^C/, Ci, T)} . (3.2) 

Here is the action for a classical path going from the point z = Q to the point z = Q 
in time T. The action functional is expected to be that appearing in the path integral 
representation of the exact propagator. The amplitude A'j.educed; the pre- exponential factor, 
is then given by a gaussian approximation to the integral over deviations from the classical 
trajectory. Such a semiclassical approximation should be accurate when j is large. 

If a continuous-time path integral is "derived" by inserting intermediate overcomplete- 
ness relations into ( p.l| ) and taking a formal limit A^ oo, then we find I^Sj 



K{Cf. a T) = dfiiz, z) exp{S(2(t), z(t))}, (3.3) 



where the path measure dn is 



^ 27 + 1 d'^z 

dij^m. m = n ^ ."^2.+2 ^ (3-4) 



and the action S{z(t), z(t)) is given by 

S{z{t), z{t)) = J {ln(l + CfZ{T)) + ln(l + z{0)Q)} + [ {^^^ " ^)} ^^-S) 

Here the classical hamiltonian, H{z, z), is related to the quantum H by 

H{z,z) = {z\H\z)/{z\z). (3.6) 

The paths z{t), z(t) obey the boundary conditions z{0) = Q, z(T) = (j, but ^(0), z{T), 
being actually z{0 + e) and z{T — e), are unconstrained, and are to be integrated over 

When we regard S as the phase-space action for a classical system |HT]], the explicit 
boundary terms, which appear naturally in the discretized path integral, serve to ensure 
that both the first-order Hamilton equations and their boundary conditions are compatible 
with the action principle. To see this, make a general variation in the trajectory, including 
variations in the endpoints. We find that 
2j£(r) , 2jz(0) ,^ 

- r h (A - ) ^ {-i^ - '§) ) <^-^' 

There are no boundary contributions proportional to 6z{0) or Sz{T) because of a cancellation 
of such terms arising from an integration by parts against those arising from the variation of 
the explicit boundary terms. Equating the variation of the action to zero therefore requires 
the classical path to obey the Hamilton equations 

^ Xl + zzfdH . .{l + zzfdH 

together with boundary conditions that fix z{Qi) = Q, and z{T) = Cf- 

The quantities z(0) and z{T) are not fixed by the boundary conditions, but can be found 
by solving the equations of motion. If we know the action for the classical path, they can 
also be read off from the Hamilton- Jacobi equations that follow from (|3.7|) , viz: 

dS,, 2jz{T) dS,i 2jz{0) 
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(3.9) 



In general ^(0) will not be the complex conjugate of z{0) = nor will z{T) be the complex 
conjugate of z{T) = C/- This means that if we write z as x + iy and z — x — iy, then, except 
in special and y are not real numbers. 

The Hamilton-Jacobi relations also tell us that 

t = t-' (3.10) 

showing that Sd is a holomorphic function of d, and an anti-holomorphic function of C/. 
These analyticity properties of Sd coincide with those of K. This is reasonable since exp 5'ci is 
the leading approximation to K, and we would expect analyticity to preserved term-by-term 
in the large j expansion. Finally 

-^^-iH{Cf.z{T)). (3.11) 

The leading semiclassical approximation is exact when the quantum hamiltonian H is 
an element of the Lie algebra of SU{2). For example, if ^ = cuJ^, then 

^ ~2'2' — 1 

H{z, z) = N\z\H\z) = ^3^^ (3.12) 

and 

dH 2jujz dH 2jujz 



dz {l + zzf' dz {l + zzf 
The equations of motion are therefore 



(3.13) 



z — icuz, z — —iujz. (3-14) 
The solutions obeying the appropriate boundary conditions are 

z{t) = e--*G, z{t) = e-(*-^)C/, (3.15) 

so 

^(r) = e— ^Ci, ^(0)=e--^C/. (3.16) 
It will only be in exceptional circumstances that z{T) — {(f)* or ^(0) = (Ci)*. 
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Inserting the solutions (|3.15| ) into the action we find 



5ei(C/, 0, T) = J {ln(l + (fQe-'^^) + ln(l + C/Oe"^-^)} + [^J^YT^ 



This is to be compared with the exact propagator 



2zz . . zz — 1] 
+ zz zz + 1 ) 

(3.17) 



K = (C/|e-^^^|0) = e^"^'^(l + = exp 



cl- 



(3.18) 



When the hamihonian is a more general element of the enveloping algebra {i.e. a poly- 
nomial in the generators) there will be corrections to this simple result. 

IV. FLUCTUATION DETERMINANT 



The prefactor in the semiclassical propagator comes from integration over gaussian fiuc- 
tuations about the classical trajectory. To evaluate these, we consider the second variation 
of the classical action, holding z{0) = Q and z{T) = C/ fixed. We will write 

1 



where 



Here 



S = S^i + 6S+ + 



6'S = -I 



2j 



il + zz) 



{5z 5z) 



-idt + A B 
B idt + A 



'5z 
6z 



dt. 



^ _ 1 ^ d {l + zzfdH ^ d {1 + zzY dH' 



B 



2\dz 2j dz dz 
(9 (1 + zzf OH 
dz 2j dz ' 
-_ d {l + zzfdH 
dz 2j dz 

When z{t), z{t) are the classical path, then 6S = 0. 



2j dz 



(4.1) 



(4.2) 



(4.3) 



On making a change of variables 
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Sz = {1 + zz)ri, 

6z= {l + zz)f], (4.4) 



we see that we have to compute the quadratic path integral 

. _ ( fTl _ \-idt + A B 
i^reduced OC / d[ri]d[ri] exp I -2ij -(v v) _ 

I Jo 2 B idt + A 



V 



dt } . (4.5) 



V 



-iasdt + M (4.6) 



This path integral is proportional to Det '^V, where the matrix differential operator 

-idt + A B ' 
B idt + A_ 

is subject to the boundary conditions //(O) = and f]{T) = 0. (We will use the symbol 
"Det " for functional determinants and "det " to denote the determinant of a finite matrix. 
Similarly "Tr" and "tr".) 

There are several subtleties involved in calculating Det T>. The most obvious is that the 
boundary conditions imposed on V are not in the class that make it self adjoint. Although 
V and are formally the same differential operator, self-adjointness requires, in addition, 
that their domains of definition coincide ||30|. It is not hard to see that the only boundary 



condition on V that leads to an identical boundary condition for is 77(0) = e^^'^f]{0) 
and ri{T) = e^^'^fi{T) for some real angles ^r- Indeed if 5 = i? = 0, then V with our 
boundary conditions has no eigenf unctions — never mind a complete set. The determinant 
cannot be expressed as an infinite product of eigenvalues, therefore. Diagonalizability is not, 
however, a fundamental requirement for defining a determinant. There exists a well-defined 
Green function G = and we should be able to obtain the determinant by varying the 
parameters and using the identity 5 In Det P = Tt {V^^6V}, which holds even if V is not 
diagonalizable. 

A potential pitfall in this approach is that the variation 6 In Det T> is given by 

SlnDetV = Ti{V-^5V} = ti {G{t,t)5M}dt, (4.7) 

Jo 

but the Green function G(t,t') is discontinuous at t = t'. We might have a different expres- 
sion for the variation depending on whether we choose to evaluate G{t, t) as G{t, t + e) or as 
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G(t, t — e). The jump in G is, however, proportional to (J3, and tr {a^^SM} = 0, so we have 
reason to hope that there is no actual ambiguity. 

If we agree to interpret Git, t) as ^{Git, t + e) + G(t, t — e)), then the formal calculation 



is straight forward ||32|, and we merely summarize the results: 
We begin by defining the matrix 

.'Vrit) Vo{t) 



(4.8) 



.VHt) 7]o{t) 

Here the column vector {r]o{t), f/o(t))"^ is a solution of T'\l/ = obeying the boundary condi- 
tion ?7o(0) = 0, "^0(0) = 1, and {ririt), fjT{t)Y is a solution with rixiT) = 1, fjTiT) = 0. The 
determinant of $(t) is an analogue of the wronskian and is independent of t. We find that 
DetP = Cdet $, where G is some constant independent of H. 

Since det $ is time independent, we may conveniently evaluate it at t = T, where 

1 VoiT) 

or at t = 0, where 

Vt{0) 
VAO) 1 

By relaxing the conditions that r/(T) = fi{0) = 1, we may interpret these results in terms 
of the variation of the endpoints of the classical trajectory as we vary the initial points. 
That is 



G-^BetV 



mm, 



(4.9) 



G-^DetV 



vAo). 



(4.10) 



G-^DetV 



\dr]{T) J \dr]{0) J 



(4.11) 



or, in terms of the original variables. 



-1 



(4.12) 



The equivalence of these two expressions for the determinant is not immediately obvious, 
but from the Hamilton- Jacobi relations 

dS,, 2jz{T) dS,, 2jz{0) 



f 



l + CfziTY dQ 
12 



1 + ^(0)0^ 



(4.13) 



and the equality of mixed partials of 5'ci, we obtain 

d^S^i _ 2j dz{T) _ 2j dzjO) 
d0ff ~ (1 + C/^(T))2 dQ ~ {l + mdr dCf ' 

Both expressions in ( |4.12| ) thus reduce to 



(4.14) 



CBet-^V ^ il^lMmi^im^. (4.15) 
Our calculation of the fluctuation determinant suggests, therefore, that 

;.-„,(C„c.,T, I ( '^^^^-'^'^f^^'°'^-' ||-)%.p5.(0,C..r). (4.16) 

(The proportionality constant is fixed by the requirement that this expression reduces to 
(C/IC) whenT = 0). 

As indicated by the "?" over the equals sign, there are problems with this expression, 
and it is not quite correct. 

The first problem is that, despite the optimism expressed above, there is a degree of 
indeterminacy in the calculation of the functional determinant. To see this, make the sub- 
stitution 

^(t) ^ e-'^Wr/(t). (4.17) 

in the path integral ( [4.51 ). The measure is unchanged, but we replace V with where T> is 
the matrix operator T) with 

A = A + dtO, 
B^B = e-2^^W5, 

B = e^'^^'^B. (4.18) 

The value of the path integral must be unaltered by this change of integration variables, but 
the solution to 
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(4.19) 



-idt + A B 1 fr,{t) 
'B idt + Al \f]{t) 

with 7^(0) = 0, r/(0) = 1 is now (e-*(^W-^(°»r/o(t), e^(^W-^(°))r7o(t))^. The determinant, as we 
have calculated it, is therefore not invariant, but ends up multiplied by e"**^^*^"^-*"^*^"^''. Our 
expression for the functional determinant has an "anomaly" therefore. 

The anomaly arises because the argument we made about the harmlessness of the dis- 
continuity in G depends on our defining G{t, t) as G{t, t ± e) with the same choice of sign in 
front of the e in both entries in the trace. If we examine the discrete version of path integral 
we see that, on the contrary, one of the entries should be evaluated with a plus, and one 
with a minus. Our calculation of the determinant assumed that we could interpret G{t, t) as 
\{G{t,t + e) + G(t,t — e)), so our formula for the determinant is only correct if both terms 
in tr {a^SM} are separately zero. This will only be the case for operators T> with A = 0. 
Fortunately the discrete path integral does permit the change of variables described above, 
and we may use this freedom to force the diagonal entries. A, to zero before computing the 
determinant. The correctly regulated functional determinant therefore differs from its naive 
value by a multiplicative factor. 

Including the correction to the fluctuation determinant, the semiclassical propagator 
becomes 

(4.20) 



where 



^^'''^-2[dz 2j dz^dz 2j dz r ^^-^^^ 



is the coefficient appearing in ( |4.3|) 

The manoeuvre of setting A to zero before evaluating the fluctuation determinant appears 
(although without explanation as to why it was necessary) in the previously cited paper by 



Kochetov that provided part of the motivation for our present work. Kochetov therefore 
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gets the corrected expression ( [4.20|) . It seems, however, that the "extra phase" (it is a phase 
only in the simplest cases), ^ Jq A{t)dt, was first obtained by Solari from a careful 
evaluation of the discrete determinant. Solari also pointed out the necessity of a similar 
correction in the harmonic oscillator coherent-state path integral, which has a fiat phase 
space. Kochetov's discovery of the correction seems to have been independent of this earlier 
work. 



Because of the extra phase, ( [4.20|) gives the correct, indeed exact, semiclassical prop- 
agator for the case H = uJ^, and also for any hamiltonian consisting of (possibly time 
dependent) elements of the Lie algebra of SU{2) [^ . 

V. SHORT TIME ACCURACY 



The Solari-Kochetov phase also solves a second problem with ( |4.16|) . In contrast to 



the configuration space propagator, which diverges as T~2, the coherent-state propagator 
K{(f, Q, T) is analytic in T near T = 0. This is because of the finite spread of the coherent- 
state wavefunctions. To first order in T we have 

K{Cf,aT) ^ (C/|e-^^^|0) ^ {Cf\Q-^T{Cf\H\Q 

= {Cf\Q{l-^TH{CfXi))■ (5.1) 

(In the last equality we have exploited analyticity to observe that the off-diagonal {Cf\H\Q), 
is obtained from the diagonal (C|-f^|C) by the the simple replacement C ^ 0? C ^ C/ •) 
Now, from the Hamilton- Jacobi equation 

dT 

we have 



-zH{Cf,z{T)), (5.2) 



5d(C/, 0, T) = 5ci(C/, 0, 0) - iTHiCf, 0) + 0(T2), (5.3) 

while 
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S^Cf, Q, 0) = 2j ln(l + CfQ) = In (C/|0). (5.4) 

Thus, in order to get agreement between ( [4.201 ) and (|5.1| ), the fluctuation determinant must 
make no 0(T) contribution to the propagator. A short calculation shows, however, that 

Fortunately this contribution is exactly cancelled by the 0{T) contribution from the Solari- 
Kochetov extra phase. 

We now ask how well does the semiclassical propagator do at next order in the short- 
time expansion. In order to provide a systematic grading for the terms, we will regard the 
hamiltonian H as being 0{j). The entire action is then homogeneous of degree one in j. 
With this assumption, and by analogy with the usual semiclassical expansion in powers of 
h, we expect that 

K{(j,Ci,T) = /^reduced " exp{5'cl} ■ 

where Sc\ is 0{j), while the prefactor, -ft'reduced, is 0{j^). 

At short time the exact coherent-state propagator is certainly of this form. To demon- 
strate this, expand 

(C/|e-^^^|0) = (C/IO) - ^T{Cf\H\Q - ^(C/li^^lO) + ■ • • • (5.7) 

Now {(f\H\Q) = {(f\Q)H{(jr,(i), but some work is needed to evaluate {Cf\H'^\(i). 
Inserting an overcompleteness integral, we have 

= ^ / (1 (1 + C/^)^^a + ^Or^HCCf, z)Hiz, Q). (5.8) 

We now perform a steepest descent expansion in the integral over the intermediate states, and 
obtain the first three terms in its asymptotic expansion in powers of 1/j. This computation 
is greatly simplified by using two shortcuts: First we need calculate only the diagonal matrix 
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(5.6) 



element {(\H'^\(). Given this, we may appeal to analyticity and obtain the general matrix 
element by setting C ^ C/ C ~^ d- Next we rotate the sphere so as to centre the 
coordinate system on the point (. Thus ( —>■ 0, and the coordinate system is locally geodetic. 
In these coordinates the saddle point of the z integral is at ^ = 2; = 0, and far fewer terms 
have to taken into consideration. 

To return to the original coordinates, we need to be able to recognize some S0{3) ~ 
SU{2) invariant combinations of derivatives and (1 + zzy factors. 
One easily establishes that, under the Mobius mapping 

a b 
c d 



az + b 

-, where 

cz + a 



e SU{2), (5.9) 



we have 



(Pz dPz' 



together with 



i + zzY n + z'z' 



, ^^^2 dfiz,z)dg{z,z) _ ^^df{z\^dg{z\z^ 
' dz dz + dz' dz' 



(5.10) 



:i+-)^^-(l+.V)^^, (5.11) 



and that the combination 



is similarly invariant. Thus, when we see the term d'^^fd^g appearing in the expansion 
about the stationary point z = 0, we realize that in the integral for the general matrix 
element (where the saddle point is at z = Q, z = (f) we should replace it by ( |5.12|) . 
Proceeding in this manner we find 



2i dCi dCf 



The three terms in braces in this expression are of 0{j'^), 0{j), and of 0(j°) respectively. 
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We may now re-exponentiate (|5.13|) as 



(C/|e-^^^ 10) = exp In (C/|0) - ^TH{Cf. 0) - t^T 



lJl + CfCi?dH OH 



2j dQ dCj 



+ ■ ■ 



X 



rp2 ^ 



4 (2j)2 V^C 
Again using the Hamilton Jacobi equation 



2 " H 



(5.14) 



dT 



tH{Cf,z{T)), 



(5.15) 



and the equation of motion for z(t), we may generate the Taylor series for Sci{T). We 
immediately verify the term in the exponential is the classical action to 0{T'^): 



5ei = ln(C/|G)-zTi7(Cy,0)-T^T 



2j ^C^ dCf 



(5.16) 



The expression in the square brackets in ( |5.14|) must be the prefactor, and is manifestly 
0(j°). It is a little tedious to verify that our formula for the pre-exponential factor, including 
the Solari-Kochetov correction, reduces to exactly this, but it is so. To collapse the terms, 
it helps to use the identity 



[l + zzf 



dzdz 



[l+^z) 



,dH OH' 



dz dz 



2(1 + ZZ? 



dH dH 



+ \{l + zz 



+ il + ZZ 



dH d 



dz dz ' \^ ' ~ ' dzdz I ' ^ ' ~' dz dz 



,1 + ZZ 



dzdz j 



which is most easily established by noting that all terms are invariant, and, at 2; = 0, both 
sides reduce to 



' 52 \ dH dH 
dzdz I dz dz 



(5.18) 



The semiclassical expression, therefore, has errors of at most 0{j ^) at short time. Our 
expectation is, of course, that it has this degree of accuracy uniformly in T. 
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VI. CONSISTENCY 



A further test of the correctness of (^4.2CI| ) is to verify its consistency under dissection of 
the classical trajectory. The exact propagator must satisfy the sewing condition 

K(Cj,Q,h + t,) = ^-l^J ^^^^^ir(C/,e,t2)i^(?,0,ti), (6.1) 

which follows from the definition of K and the overcompleteness condition ( p^.4| ). The semi- 
classical approximation to K should obey a similar condition, but with the exact integration 
over the intermediate states replaced by a suitable stationary phase approximation. 

Since Kgd ~ exp Sd, we begin with the relationship between the action for the total path 
from (i to (f, and the actions for the two segments from Q to the intermediate point ^, and 
from ^ to (f. To eliminate the redundant intermediate-point boundary terms we must define 

SiCj, Q, h + h) = S(Cj, e, h) + S{1 0, ti) - 2j ln(l + ee). (6.2) 
We will write this compactly as 

S'tot = ^2 + ^i-2jln(l + eO- (6.3) 



In writing ( |6.2| ) we have tacitly assumed that our chosen starting ^ of the second path 
segment coincides with the dynamically determined endpoint z{ti) of the first path segment, 
and that the dynamically determined starting z{ti) of the second path segment coincides 
with our chosen ^ endpoint of the first path segment. This will not generally be the case — 
but it will be when ^, ^ obey the stationary-phase equations 

^ = ^ = 0. 6.4 

Taking into account the analyticity properties of 5*1 and 5*2, these are 







1 + er 



di i + a ^ ' 
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Comparing ( |6.5|) with the Hamilton- Jacobi equations confirms that = -2(^1) and = ^(ti), 
where ^c, C is the stationary phase point. 

To evaluate the integral over small deviations from the classical stationary phase point, 
we set ^ = + ^, ^ = C + V- We expand 



'S'tot — '5'toti 



2j 



2 (1 + UcY 



1 —a 
-13 1 



V 



(6.6) 



where 



a 



d 



2J «r«="2j(^"«='c'^"«='C 



(6.7) 



and 



(6.8) 



(The second equality in these equations uses the stationary phase equations.) 

We now put together two semiclassical propagators and perform the gaussian integral 
over the deviation from the stationary phase point. Using the semiclassical Solari-Kochetov 
form ( |4.2CI| ) for the propagators on the right-hand side of Eq. ( |6.1| ), we get (with T = ti -1-^2), 



K, 



2j + l 



comb 



■ exp iSi + S2- 2j ln(l + ^e^c) + / Adt- -S'S 



X 



2j 



dCfd^c 



2j 



(6.9) 



Notice that, as with consistency test of the ordinary Feynman path integral |^3[, the measure 
and the prefactors, including the Solari-Kochetov "extra-phase" term, are all being treated 
as constants. The integration involves only the variation of the classical action 



2j 



iv, v) 



1 —a 
-(3 1 



(6.10) 



and yields, along with other factors, the inverse square-root of the determinant 

1 —a 

D = 

We now use the result, established in the appendix, that 



(6.11) 
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dCfdQ 2j dCfdi,dl,dQ 

to obtain 



1 —a ^ 



-(3 1 



(6.12) 



(6.13) 

The semiclassical approximation therefore reproduces itself except for a nigghng factor of 
(2j + l)/2j, which is due to a conflict between the normalization of the measure and the 2j 
appearing in the exponent. 

Although this discrepant factor approaches unity in the large- j limit, it is nonetheless 
disturbing. Each of the infinitely many gaussian integrations that constitute the semiclassical 
approximation to the path integral ought to be indistinguishable from our single gaussian 
integration over the intermediate point ^. We should, therefore, be able to dissect the path 
into arbitrarily many parts without affecting the final answer. This is not currently so, and, 
in particular, the limit of large j does not commute with the limit of a large number of 
intermediate points. 

The origin of the discrepancy is not hard to find. In the large-j limit the effective radius 
of our spherical phase space becomes large, and, near z = 0, the spin-j reproducing-kernel 
relation 



2j + 1 f d^z 



+ = (C2IC1), (6.14) 



TT J {1 + zzy 

or more explicitly, 

^ / 7T^ii + zzr\i + C2zr\i + zCiY' = (1 + C2Cl)^^ (6.15) 

TT J [i + ZZ) 

should contract to a suitably scaled version of its flat-phase-space analogue 

— e-^"e^2^e^^^ = e^^fi. (6.16) 

71 



Because it is a gaussian integral, the leading stationary phase "approximation" to ( |6.16D is 
exact. 

21 



If we make the obvious large j estimates 

(1 + zz)-^' ~ e-'^~^ (1 + C^zf^ ~ e2^'^2^ (1 + zO)2i ^ e^^^^\ (6.17) 

while regarding the sphere measure (1 + Iz)'"^ as a prefactor, we do not get exactly 

2j 



TT 



but instead (2j + l)/2j times this. 

If we keep terms higher order in l/2j, both those coming from the measure and those 
from going beyond the quadratic approximation to the exponent, they will of course correct 
the error. What we really need, however, is a partitioning of the integral on the LHS of 
( |6.15| ) such that the leading steepest descent approximation will agree with the RHS. This 
will happen if regard the expansion parameter as 2j + 1 and not 2j. To see this, break up 

/ = ^ / jT^.a+z^r'^ii+c^zr^ii+Kir^ (6.19) 

as 



/ = ^ / jr^,9~\z,z)e^'^^'^'^^^^^-^ (6.20) 

TT J [i + zzy 

with 

g{z, z) = {l + zz)-\l + C2^)(l + ^Ci)- (6.21) 
The critical point of the function in the exponential is at ^ = (^2; ^ = Ci? ^ind 

^?(C2,Ci) = (l + C2Ci), (6.22) 
d'^lwg 1 



dzdz 



(6.23) 



Thus 



^ (l+C2Cl)^ 
= ^(l + C2Cl)^-^-y^(l+C2Cl)^ 

= (l + C2Ci)'^' (6.24) 
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The leading term of the asymptotic expansion of / in powers of l/(2j + 1) is therefore exact. 
This observation suggests rewriting the semiclassical approximation to K as 



i^scl(C/,0,T) 



v/2jTT \dCfdQ 

where Sc\ = (2j + l)S'ci/(2j), and 

ifil + zzYd^H 



^ explSci{Cf,Ci,T) + ^ Qdt 



(6.25) 



+ H(z, z) 



(6.26) 



j \ 2 dzdz 
is the term required to make (|6.25|) numerically equal to (^4.20|) . 

With this repartitioning of terms between the exponent and the prefactor we have exactly 
the same classical equations of motion, but now 



comb 



:i + C^c) exp ^ 5i + ^2 - (2j) ln(l + C^c) - \^''~S 



X ■ 



(2j) Uc/^^cC^o 



exp 



2 7o 



Qdt 



(6.27) 



where 



(6.28) 



and 



2i 



1 

-/3 



—a 
1 



(6.29) 



and (|6.8|) by putting tildes on S2, 



The quantities a and /? are obtained from Eqs. 
and j. Note, though, that a = a, and P = P- Note also, that we have inserted a factor of 
(1 + C,c^c) in the integral to compensate for the extra factor of (1 + that was taken from 
the measure into the exponential to complete Stot- Thus part of both the measure and the 
prefactor are varied in determining the stationary phase, and get integrated over, while part 
is regarded as a constant. 

The integration in Eq. ( |6.27] ) can be done at once by noting that all equations in Appendix 
A are unchanged if we put tildes on the actions, j, a, and f3 everywhere. In particular, the 
identity ( [A^ ) holds with tildes. We thus obtain 
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1 

all unwanted factors of 2j + 1 and (1 + C,c^c), having cancelled. Thus, with this form of 
stationary-phase integration, the propagator reproduces itself exactly. 

What this means is that the semiclassical approximation must be tacitly using ( |6.24| ) 
in making each of the many integrations that go into the gaussian approximation to the 
path integral. Once we realize this, we see that there is no need for the mysterious diver- 
gent normalization factor, Af = limAr^oo(l + that plagues most treatments of the 
semiclassical spin propagator. 

The appearance of j + 1/2 as the large parameter in the fluctuation integral has been 
remarked on before by Ercolessi et al. 0] and by Funahashi et al. |^5[. The former worry 



that it is inconsistent to include fluctuations of the measure in the gaussian integral without 
also considering their effect in the saddle point equations. In our case all terms that are 
being integrated over do appear also in the equations determining the saddle point. 

Note that the correction Q vanishes for Larmor precession where H = ujj-^. In this case, 
as we have seen earlier, 

5-^1 = 2j ln(l + C/Oe-*"^) + ijujT. (6.31) 
5* is obtained from this by the substitution j ^ j + |, so 

, (6-32) 



d(fdQ (i + oc.e--^)2- 



Thus 

1 / d''S 



1 

2 c \ 2 



= e*"^(l + C/Oe-'"'^)'^ (6.33) 



which is the exact answer. 
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VII. AN EXAMPLE: H = uJ^ 

As an application of the semiclassical formalism consider H = uJ^. This hamiltonian 
is time reversal invariant, and we might worry that a hidden shift j — >■ j + 1/2 would 
compromise the Kramers degeneracy expected when j is half integral. 

The classical hamiltonian corresponding to H = z/J| is 

This should be compared with the "naive" classical hamiltonian 

Hn.^.e = (^^)\ (7.2) 

\zz + 1/ 

which is what we would get if we simply expressed the classical direction-dependent energy 
z/j^ cos^ 6 in terms of the stereographic coordinates on S*^. 

The hamiltonian (|7.1| ) leads to the classical equations of motion 

z = iuj(z,z)z, z = —iuiz, z)z, (7.3) 

where, with fj, = uj{j — 1/2), 



j J \zz + 1. 

Since these equations imply the time independence of the product Iz, uj is itself time inde- 
pendent and the solutions may be written down directly as 

z{t) = e-^-*0, z{t) = e-(*-^)C/- (7.5) 
Here uj is to be determined by the self-consistency condition 

As we will see below, this equation has an infinite family of solutions. Here, we wish to 
consider how various quantities scale with j. By demanding that Eqs. ( |7.3| ) continue to be 
meaningful as j — 00, we see that we must have /i = 0{j), u = 0(1), and u = 0(1/ j). 
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The classical action for the solution ('[7.51) is 



5ei(C/,0,T) = 2jln(l + e-^-%0: 



2j ln(l + e-^-^CfO) + ^T{jcu + ^uo^ - Uv}. 



(7.7) 



4/i 2- 

The apparently cosmetic rewrite in the last line leads to a useful way of looking at the 
problem. Define 



S^iCf, 0, T) = 2j ln(l + e-'^^CfQ) + iT{ju + ^u'} 



7 2 

J , .21 



4/i 



(7.8) 



where we regard uo as an independent variable. The equation 



- = iT] 



e , ( j 



(7.9) 



then shows that the consistency condition on u is equivalent to dS^/du = 0. We can also 
use Si^{( to express the second variation of Sd required for the prefactor A. By 



differentiating the Jacobi equation ( p.9|) we have 



2j dz{T) 



dCfdQ (1 + C/^(T))2 dQ 

and from this we find, with Eq. ( [7.5|) , that 



(7.10) 



g^gel(C/,C»T) 



2i 



We now differentiate the condition dS^^/du = with respect to Q. This yields 



(7.11) 



+ 



0. 



Using this result to eliminate (du/dQ) in Eq. ( [T.lll ), we find, after a little algebra, that 



(7.12) 



;i + C.^(T))2 2/i 



(7.13) 



Substituting Eqs. (|7.5|), (|7.7|), and ( |7.13| ) into the basic semiclassical form ( [4. 201 ) for the 
propagator, we obtain 



26 



iTf-V'^ /S=S„\"'" f„ iT, . , i 



A- = i:h7 exp S„--(. + ..)-H-/^ . (7.14) 



T 



The sum over u is to be performed over all solutions to Eq. ( |7.6|) . 

The utility of Si^{(j,(iyT) is not hard to understand. We are trying to evaluate 



while we already know that 



rn^-j [j + m)\{j - m)\ 



m=-j 



expSu;oi(f,(i,T), (7.16) 



where 



^-o(C/,0,^) = 2jln{l + e-'^^CjQ)+zTju. (7.17) 
From the identity 



V A-nu J-oo ^ 
we have the exact relation 



■ TT I T f \ — uo'^T 



4z/ 



V4^ / ^^exp|2jln(l + e-^-%0)+^nJ^ + ^}|- (7-19) 



Given the form of the classical action ( |7.7| ), that /i ^ j z/, and the occurrence of 



{d'^Sijj/ duj'^)~^f'^ in the prefactor, it is clear that the semiclassical approximation is attempt- 
ing a stationary phase approximation to this integral over uj. That this approximation is 
indeed indicated can be seen by evaluating (d'^S^/duj^). From Eqs. ( [7.9| ) and ( [7.6|) , we find 
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which scales as j as j — oo. 

We now write the exponent in Eq. ( |7.19| ) as S^^ — iTjuj'^/8fi. Since the second term is 



0(j°) as j —>■ oo, we may regard it as part of the pre-exponential factor in carrying out the 
stationary phase integral. In this way, we obtain 

The pre-exponential factors in the preceding equation agree with those in Eq. (|7.14|) to 
terms of order unity. To see whether the exponents agree, we must discuss the effect of the 
Solari-Kochetov phase. We find that 

f d {l + zzf OH d{l + zzydH\ 

~ 2 vai 2j 'dz ^ 2j aFy 

4/i zz 



J (i + zzy 

,2 



The term in parentheses serves to cancel [up to 0(1)] the second term in the exponent in 
Eq. ( |7. 14| ) , and the ju'^/ifi term serves to correct Suj as needed in Eq. ( [7.21| ). Thus our 



semi classical formula is indeed accurate up to 0(1) as j —>■ oo, and we may be confident 
that spectral properties (Kramers degeneracy in particular) derived from it by constructing, 
say, the Green's function or density of states, will be faithfully given. 

Having demonstrated the formal equivalence of i^sci and inexact, we turn to the actual 
nature of the solution. Let us first rewrite the self-consistency condition (|7.6| ) as 



2iCoT + In ( ) = In a, (7.23) 

I — U) 



where u = juj/2fi, r = fiT/ j, and a = (fQ. In the limit r ^ oo, the left hand side of (|7.23|) 



must remain finite, suggesting that ~ 1/r. A development in powers of 1/r shows that 
we may write 

Cu^ ^ / +0r-^. 7.24 
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Since no restriction has been placed on which branch of In a is to be taken, this solution is 
infinitely multivalued, as asserted above. To leading order in 1/r, different solutions differ 
by additive amounts nvr/r, where n is an integer. 

On the other hand, at r = 0, Eq. ( [7.23| ) has a unique solution, u) = (a — 1)/(q; + 1). 



The apparent contradiction with the earlier argument for an infinite number of solutions is 
resolved by noting that if, as r 0, we allow to diverge as 1/r, the left hand side of 
( [r.23| ) again remains finite. Another development in powers of r reveals that 

% ln(— a) 2 8i 
2 r ln(— a) [ln(— ajj'^ 

which is also multivalued on account of the infinitely many branches of ln(— a). 

We can gain further insight into the nature of the propagator and the values of lo at 
the relevant stationary-phase points by working with initial and final states on the equator 
of the sphere: = c**^', C/ = e~^'^f . When j is large, the problem should be essentially 
equivalent to a massive particle constrained to move on a ring of circumference 27r. If we 
write the hamiltonian for the latter as L^/2M, where L is the orbital angular momentum, 
and M the mass, we expect the results for the two problems to be similar with M = 2v. 

We start by considering the propagator for Larmor precession. Employing the leading 
large-j estimate 

2j! 22^- 



e"™ l\ (7.26) 



{j + m)\{j - m)\ 0rj 
and using the shorthand A0 = (pf — (pi, we may write 

^tlj (j+m)!(j-m)! 

If T ^ j^^/^/u;, the summand will have widely varying phases over the range of m values 
that contributes to the sum, |m| ~ y^. By extending the sum over m to infinity and using 
the Poisson summation formula (taking care that m takes half-integer values when j is half 
integral), we find 
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n 

where the (—1)" factor is present only when j is half integral. This form is better suited to 
studying the large j limit (for fixed T). In that case, (|7.28|) , regarded as a function of tu, is 



sharply peaked aX uo = uOn = (27m — A(j))/T. These are the angular frequencies that allow 



uniform precession between 0j and in time T. We now recall that Eq. ( 7.28 ) is nothing 



but exp(S'(^o)- If we substitute this form into Eq. (|7.19|) , and take into account the factor 
exp{iu;^T/4z/} in determing the saddle-point frequencies, we find that they become complex 

UJri — UJn I 1 



uTj 



OOn + (7.29) 



Not surprisingly, this is just what we found in Eq. ( |7.24] ). The result reflects the fact that, to 



move at the required speed, the hamiltonian trajectories must move off the equator. There is 
then no real trajectory between the classical endpoints, and we must exploit the freedom to 
have trajectories where C/ 7^ z{T)*. When j is large, however, Hamilton's equations provide 
large velocities dose to the equator, and the imaginary parts of uj are correspondingly small. 



Performing the integration in Eq. ( [7.19|) , we find 



^g«0/|g-*fTJ?|g«0>^ ~ o'i.i^-i.iM ^ ^ ^-|(A<ji-27rn)V(l+uVT) ^ (_1)"^ (7.30) 

{l + ijuT)2 n 

where, again, the last factor is only present when j is half-integral. This form should be 
compared with that for the massive particle p6| 



We have incorporated an Aharonov-Bohm phase $ into the result. This phase should be 
vr when we compare with half-integer spins, and the resulting pairwise degeneracy of the 
energy levels is the particle-on-a-ring analogue of Kramers degeneracy. 

The similarity between Eqs. ( |7.30|) and (|7.31| ) is evident. Notice how j sets the time 



scale for the crossover between the large-T regime, where the spin behaves essentially as 
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a particle of mass 2v on the ring, and the short-time regime where the finite range of the 
coherent-state wavefunctions cuts off the 1/\/T divergence. 

Note that we have ignored the difference between /i/j^ and v in the above comparison, 
since as discussed while showing the equivalence of -K^sci and -ft'exact, the error incurred is 
of order relative to the leading term in the action. The semiclassical approximation 
therefore correctly obtains the first two terms in the large-j expansion. 



VIII. DISCUSSION 

In the previous sections we have used the continuous-time path integral to motivate 
a semiclassical approximation to the coherent-state propagator for spin j. Although our 
derivation of the semiclassical propagator is purely formal, and the resulting expression 
must initially have only the status of a conjecture, we have demonstrated its correctness by 
verifying its short-time accuracy to O(T^), and checking its consistency under dissection of 
the path. From these two properties we may conclude that our expression is accurate to 
0(j*') uniformly in time. 

In our derivation it was necessary to take into account an "anomaly" in the evaluation 
of the functional determinant of the Jacobi operator. This is the only place where we had 
to appeal to details of the discrete version of the path integral. Regulating the determinant 
in a manner consistent with the discrete path integral results in a correction to the naive 
expression for the prefactor. This correction had been noted before, by Solari p3 and by 



Kochetov ||2^, but its importance does not seem to have been widely appreciated. 

We have also discussed an example where an infinite number of classical trajectories 
contribute to the propagator. Here we again saw how the Solari-Kochetov factor is essential 
in obtaining the correct result. 

A calculation of the Solari-Kochetov correction to the tunnel splitting between classically 
degenerate spin states will be reported in a separate publication. 
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APPENDIX A: COMPOSITION OF PATH-DENSITY FACTORS 



In the this appendix we derive ( |6.12|) . We begin by restating the stationary phase 
conditions ( |6.5D : 

dS2(Cf,o 








5j 



2je 



Consider how the first of these evolves as we vary (j-. We find that 
d fdS2 2jC \ 



(Al) 







d^S2 



2jC d^c 



2j d^. 



dCfdi, dec (1 + ec^c)^ (1 + ac)2 

d'S2 , die I d'S2 , 2jll \ 2j C 



+ 



+ 



(A2) 



dCfd^c dCf \ oe (1 + Uc)'J (1 + Uc)' dCf 
In the last line, we recognize the expression in parentheses to be 2j/?/((l + ^^icY, where 
j3 is the coefficient appearing in (|6.6| ). By differentiating each of the two stationary phase 
conditions with respect to C/ and 0, we get a total of four such equations. These may be 
summarized as 

- - 2 / ^ 



1 —a 

Taking determinants, we obtain 
1 —a 



die aic_ 

da dCf 



1 / \ ^ ^ 



2j 







(A3) 



dCfd^c 



-f3 



die dic_ 

dQ acj 



(2jT d^dQdCfd^c 



(A4) 
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We now recall that the gaussian integration in Eq. ( |6.9|) leads to the inverse- square root 
of the precisely the first determinant in Eq. ( |A4| ) . This equation expresses this determinant 
in terms of the second derivatives of Si and 52, and the jacobian d{^c,^c)/d{Ci,Cf)- The 
derivatives of 5*1 and 5*2 will cancel with the prefactors in Eq. ( |6.9| ), leaving only the jacobian. 
We therefore turn to its evaluation, and show that it can be written in terms of the second 
derivatives of ^tot with respect to (j and Q. We express iStot as 

5tot = 52 + 5i-2jln(l+^0, (A5) 
and take note of the fact that both C,c and vary as we vary (j- and Q. We have 

dCfdQ \ di, dQ 1 + dQ 9C dCi dQ 1 + dQ J 



dCf \dQ {die 1 + a J V \di, 1 + J 

The expressions in braces in the last line are the stationary phase conditions, so they are 
zero, as are their derivatives. Thus: 

S^^tot d (OS A d^Si C 



(A7) 



Taking note of the fact that the derivative of Si with respect to Q is at fixed while we 
have useful expressions for the derivative including the variation of we interchange the 
order of differentiation, and write 

eyiot ^dL{d_ (dsA _ d^Si c 

C d ( 2ji, \ dl.d^Sidl, 



dCfdCi \l + Uc) dCf dt 9Ci 
C d ( 23i, \ dl, d ( 23i, 



In going from the second line to the third, we used one of the equations from 
Putting this together with ( [A4|) yields 



(A8) 
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which is identical to Eq. 



1 



-a 



(A9) 
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